The peculiar velocity field: constraining the tilt of the Universe 
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A large bulk flow, which is in tension with the Lambda Cold Dark Matter (ACDM) cosmological 
model, has been observed. In this paper, we provide a physically plausible explanation of this bulk 
flow, based on the assumption that some fraction of the observed dipole in the cosmic microwave 
background is due to an intrinsic fluctuation, so that the subtraction of the observed dipole leads to 
a mismatch between the cosmic microwave background (CMB) defined rest frame and the matter 
rest frame. We investigate a model that takes into account the relative velocity (hereafter the tilted 
velocity) between the two frames, and develop a Bayesian statistic to explore the likelihood of this 
tilted velocity. 

By studying various independent peculiar velocity catalogs, we find that: (1) the magnitude of the 
tilted velocity u is around 400 km/s, and its direction is close to what is found from previous bulk 
flow analyses; for most catalogs analysed, u = is excluded at about the 2.5o level; (2) constraints 
on the magnitude of the tilted velocity can result in constraints on the duration of inflation, due 
to the fact that inflation can neither be too long (no dipole effect) nor too short (very large dipole 
effect); (3) Under the assumption of a super-horizon isocurvature fluctuation, the constraints on the 
tilted velocity require that inflation lasts at least 6 e-folds longer (at the 95% confidence interval) 
than that required to solve the horizon problem. This opens a new window for testing inflation and 
models of the early Universe from observations of large scale structure. 



I. INTRODUCTION 

The bulk flow, i.e. the streaming motion of galaxies or 
clusters, is a sensitive probe of the density fluctuation on 
very large scales. Recently there have been several obser- 
vations of a large amplitude of the bulk flow on hundred 
Mpc scales, which are in conflict with the predictions of 
the ACDM model 0,0. The bulk flow is measured with 
respect a frame in which the CMB temperature dipole 
vanishes. We define this to the CMB rest frame. It is 
usually assumed that the CMB rest frame coincides with 
the matter rest frame which we define to be the frame in 
which the velocities of matter in our horizon volume are 
isotropic. 

It is possible that the two reference frames actually 
do not coincide with each other, resulting in a "tilted 
universe" [3-5]. If the inflationary epoch lasted just a 
little more than the 60 or so e-folds needed to solve the 
horizon problem, the observable "remnants" of the pre- 
inflationary Universe may still exist on very large scales 
of the CMB. As a result, there could be some fraction of 
the CMB dipole due to the intrinsic fluctuations rather 
than observer's kinetic motion. Therefore, when the ob- 
served dipole is subtracted from the galaxy peculiar ve- 
locity data, the subtraction induces a mismatch between 



the CMB rest frame and matter rest frame. We define 
the intrinsic CMB dipole to be the CMB dipole mea- 
sured in the matter rest frame. The recent finding of 
such a mismatch between the directions of the observed 
CMB dipole and reconstructed velocity dipole [H, 0] sug- 
gested the possible existence of the intrinsic CMB dipole 
on the sky. 

Pre-inflationary fluctuations in a scalar field may pro- 
duce an intrinsic dipole anisotropy. In this scenario, the 
observed CMB dipole would be a sum of the dipole from 
our motion in the matter rest frame and the intrinsic 
CMB dipole caused by a large scale perturbation. This 
intrinsic dipole can be produced by a large scale isocurva- 
ture perturbation , but not a large scale adiabatic per- 
turbation Q . Note that our local motion relative to the 
matter rest frame is caused by small scale inhomogeneity 
(up to about the 100 Mpc scale) and will be negligibly 
affected by the very large scale (3> 10 Gpc) perturbation 
that would cause a tilt effect. 

In this paper, we estimate the tilted velocity (u) be- 
tween the two rest frames using galaxy peculiar velocity 
data. This opens a new window on testing early- Universe 
models from observations of large scale structure. 



'Electronic address: yzm20@cam.ac.uk 



1 In our context, an isocurvature perturbation is distinguished 
from an adiabatic perturbation in that the ratios of the num- 
ber of photons to baryons and cold dark matter particles is not 
spatially invariant. 
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II. LIKELIHOOD AND MOCK CATALOGS 

In order to use the galaxy peculiar velocity catalogues 
(see the following descriptions of the data), we need to 
model the velocity of galaxies in different rest frames. 
For each galaxy velocity survey, we first subtract off our 
local motion with respect to the CMB as estimated from 
the CMB dipolc. However, if there is a non-negligible 
intrinsic CMB dipole, the CMB defined rest frame will 
not correspond to the matter rest frame and thus there 
will be a residual dipole in the galaxy peculiar velocity 
survey. To test this we estimate the line-of-sight velocity 
S n of each galaxy n in the CMB rest frame with measure- 
ment noise a n . Suppose the CMB rest frame has a tilt 
velocity u with respect to the matter rest frame, then the 
line-of-sight velocity of each galaxy with respect to the 
matter rest frame becomes p n (u) = S n — r n ^Ui, where 
f n ^Ui is the projected component of the 3-D Cartesian 
coordinate u onto the line-of-sight direction of galaxy 
n. After subtracting out the "tilted velocity" , we model 
the galaxy line-of-sight velocity with respect to the mat- 
ter rest frame as p n = v n + S n , where v n is the galaxy 
line-of-sight velocity in the matter rest frame, and 5 n is 
a superimposed Gaussian random motion with variance 
a n + °* i where cr» accounts for the 1-D small scale non- 
linear velocity. It can also compensate for an incorrect 
estimation measurement noise a n [H, [2]. Therefore, the 
covariance matrix of p n ( u ) becomes 



G, 



(v m v n ) + 6 mn (al + al) 

<(?™-v(rJ) (f n -v(r„))> + S mn (al + a, 2 ),(l) 



in which the cosmic variance term is Q 
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where 



fmn (J$) 



J (frrk) xexp (ikk- (r m - r n ) 



(3) 

which can be calculated analytically (Appendix [5J . 

Therefore, the likelihood of the tilted vector and the 
small scale velocity dispersion cr* can be written as 



£(u,o-*) = — — ; -rcxp ( -ip m (u)G m ^ n (u) 

(detG mn ) 2 V 2 



(4) 

where we fix the cosmological parameters at the WMAP 
7-year best-fit values (fi b = 0.0449, O c = 0.222, h = 0.71, 
and cr 8 = 0.801 [uj). 

We parameterize the velocity as u = {u,cos(9),4>} 
where in Galactic coordinates <fi = I and 9 = ir/2 — b. 
In order for our marginalized prior on u to be uniform 
we set Prior(u, ct») oc 1/u 2 . Then the posterior distri- 
bution of the parameters (u, c*) given the data (D) is 




FIG. 1: The SN data plotted in Galactic coordinates.. The 
red points are moving away from us and the blue ones are 
moving towards us. The size of the points is proportional to 
the magnitude of the line-of-sight peculiar velocity. "X" is 
our estimate of the direction of the tilted velocity estimated 
from the SN data. 



Pr(u,cr*|D) oc Prior(u, cr*)L(u, cr*) . Before we perform 
the likelihood analysis for the real data, we have tested 
this likelihood from mock catalogs. We input a set of 
fiducial values of (ct*,u) and do 300 simulations of the 
data. The average likelihood of these simulations for 
each parameter exactly peaks at the input values, with 
the appropriate width determined by cosmic variance, 
instrumental noise and small scale velocity dispersion. 



III. DATA ANALYSIS 

We studied five different catalogs of galaxy peculiar 
velocities to constrain the tilted velocity u and the ve- 
locity dispersion <r*. The five data sets are (see [l[ for 
the procedure of excluding outliers): 

• ENEAR is a survey of Fundamental Plane(FP) 
distances to nearby early-type galaxies [H]]. After 
the exclusion of 4 outliers, there are distances to 
698 field galaxies or groups. For single galaxies, 
the typical distance error is 20%. The characteristic 
depth of the sample is 29 /i _1 Mpc. 

• SN are 103 Type la supernovae distances [12], lim- 
ited to a distance of <150 /i _1 Mpc. SN distances 
are typically precise to 8%. The characteristic 
depth is 32 fc^Mpc. See Fig. [H 

• SFI++ [13], based on the Tully-Fisher (TF) [H] 
relation, is the largest and densest peculiar veloc- 
ity survey considered here. After rejection of 38 
(1.4%) field and 10 (1.3%) group outliers, our sam- 
ple consist of 2720 field galaxies and 736 groups, 
so we divide it into two sub-samples, field samples 
SFI++F, and group samples SFI++G- The char- 
acteristic depth is 34 h~ 1 Mpc. 

• SMAC [3] is an all-sky Fundamental Plane (FP) 
[29l ] survey of 56 clusters, with characteristic depth 
65 ft- x Mpc. 



3 



• COMPOSITE is the combined catalogs (4536 data 
in total, compiled in Q],Q) which has the character- 
istic depth of 33 ft, _1 Mpc. It is a combination of SN, 
SFI++ SMAC, ENEAR and also the samples from 
SBF [l5j (69 field and 23 group galaxies, character- 
istic depth 17 /i _1 Mpc), EFAR [lf| (85 clusters and 
groups, characteristic depth 93 /i _1 Mpc), SC [I?} 
(TF-based survey of spiral galaxies in 70 clusters, 
characteristic depth 57 /i _1 Mpc), and Willick [l8[ 
(Tully-Fisher based survey of 15 clusters, charac- 
teristic depth 111 ft _1 Mpc). 

In Fig. [21 we show the marginalized 1-D posterior 
probability distribution functions of the small scale ve- 
locity and intrinsic dispersion c, , and the velocity vector 
u. The best fit and marginalized error bars are listed 
in Table HI In panel (a) of Fig. [2J we see that different 
surveys prefer different er». For the SN catalog, there 
is a smaller er», because the supernovae light curves can 
be used to calibrate the distance measurement very pre- 
cisely, so the scatter of the line-of-sight velocity is small. 
In the COMPOSITE catalog, er» ~ 450 km/s reflects the 
average value of <7* of all of the catalogs, which can be 
treated as the average value of the small scale velocity 
and the intrinsic dispersion. Using a single <7* should 
not have a significant effect, in general, because the cata- 
logs are of typical depth > 40 /i _1 Mpc or in better units 
> 4000 km/s, and the distance indicators have typically 
20% error in the measurement which means that typi- 
cally <r n > 800 km/s. So whether <7* is 400 km/s or 600 
km/s is not of much consequence if a n > a*. 

In panel (b) of Fig. [2J we show the marginalized 1-D 
distribution of the magnitude of the u. Also, from Ta- 
ble U we can see that the SFI++ and ENEAR catalogs 
can provide fairly tight upper bounds on u, and ENEAR 
peaks at zero while the SFI++ catalogs contain zero ve- 
locity within 2<T. It should be noticed that the SN and 
the larger COMPOSITE catalog can provide a non-zero 
2cr lower bound on u, in which the latter provides the 
tightest constraints. In panels (c) and (d) of Fig. [21 we 
show the marginalized 1-D probability distribution of the 
direction of the relative velocity (cos(#), <fi). We see that 
the distribution of cos(#) and <fi are very close to Gaus- 
sian, and the four catalogs predict a very similar direction 
of the velocity (Fig. [3] and Table HI. We should also no- 
tice that the various catalogs give consistent constraints 
on the tilted velocity. 

In panel (a) of Fig. we plot 2-D contours of er» 
and it, and we find that, in the posterior distribution, er„ 
and u are not correlated. The reason is easy to under- 
stand: the "tilted Universe" velocity u describes super- 
horizon dipole modulation of CMB photons, whereas c* 
describes the sub-horizon modes for small scale veloci- 
ties and the intrinsic dispersion of each individual galaxy, 
so they come from completely different origins therefore 
without much correlation. Thus, if the instrumental noise 
a n has been underestimated, the distribution of <7* will 
be shifted to smaller values, but that will not change the 
constraints on u significantly. Thus, our constraint on 



u is pretty robust with respect to the estimate of small 
scale instrumental noise. In panel (b) of Fig. [3J we plot 
the la contour of Cos((9) and <fi. We can see that the di- 
rections of tilted velocity found by different surveys are 
very consistent with each other, and also very consistent 
with the results from Watkins et al. [l[ . 



IV. DISCUSSION 

We should notice that.the direction of the large bulk 
flow velocity found by [l[ is within the ler confidence 
level of the tilted velocity here, which therefore can pro- 
vide a physical origin of the large bulk flow in [l|. In Q 
they evaluated that the probability of getting a higher 
bulk flow within ACDM with WMAP7 parameter val- 
ues to be less than 2%. In our analysis, the bulk flow 
of a galaxy survey should be accounted for in the error 
bars by the cosmic variance term (Eq. [2]). So an alterna- 
tive explanation might be that there is a feature in the 
matter power spectrum which would increase the cosmic 
variance. However, as the shear and octupole moments 
of the peculiar velocity field are not anomalously large, 
this is disfavored by the data (l9l |. 

In addition, the direction we find is consistent with 
that found in [20] which used the dipole of the kinetic 
Sunyaev-Zeldovich (kSZ) measurements to estimate the 
bulk flow on Gpc scales, but our magnitude is lower than 
what they found. However, there is an additional level of 
uncertainty in converting the kSZ dipole into a bulk flow 
which makes it difficult to estimate the magnitude of the 
bulk flow from the kSZ dipole [2(j- In the bulk flow ap- 
proach, Ref. [2| suggests that the bulk flow comes from 
scales > 300 /i _1 Mpc. If [2(| proves correct it will be 
strong support for tilt explanation since the tilted veloc- 
ity should be the same regardless of the scale probed to 
measure it. A tilted Universe scenario was also proposed 
in [j| to explain the kSZ measurement. 

An intrinsic dipole on the CMB sky caused by a tilted 
Universe may be explained by a pre-inflationary relic 
isocurvature inhomogeneity. If inflation lasts for only 
a few e-folds longer than required to solve the horizon 
problem, the scales that were super-horizon at the ini- 
tial point of inflation are not very far outside our current 
horizon today. Inflation requires that there was initially 
a fairly smooth region of order of the inflationary Hub- 
ble horizon. A physical scale of such a region at the 
onset of inflation I = ePH~ x will have the physical scale 



L 



3 w- 1 



today, where P = p + N - N min (N n 



the minimal number of e-folds to solve the horizon prob 



In 



lem and is generally assumed to be around 60) 
the short inflationary case where N is not much larger 
than TVmin, a remnant of a pre-inflationary Universe still 
exists on super-horizon scales which may be detectable. 
The quadrupolc effect, aka the Grishchuk-Zel'dovich ef- 
fect [21j , is also part of the effect of this large scale inho- 
mogeneity. 

An isocurvature perturbation can produce the "tilted 
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FIG. 3: (a) The 68%, 95% and 99.7% Bayesian confidence interval contours for parameters a*-u. (b) The 68% contours for 
Cos(8)-<fi. The black dot is the direction of bulk flow found by [lj. 
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Catalogs 


Characteristic Depth (h *Mpc) 


a, [lOOkm/s] 


u [lOOkm/s] 


£ (degrees) 


b (degrees) 


A TV {2a) 


ENEAR 


29 


2.8 ± 0.3 


q+2.2 


987 9+68.9 


n Q + 37.5 
°-°-36.0 


AN > 7 


SN 


32 


1.8 ±0.4 


4 -°-1.9 


284.91^-9 


1 q+18.8 


6 < AN < 9 


SFI++G 


34 


5.iolS:5 


n 1+1.6 
°- x -1.9 


289.7 ±34.9 




AiV > 6 


SFI++ F 


34 


6.6i°; 3 2 


n q+1.2 
z '°-1.6 


276.8^39.0 




AiV > 6 


SMAC 


65 




5.9 ±2.4 


263.8l?fj 


1 1+18.6 

^-le.o 


6 < AN < 8 


COMPOSITE 


33 


4.81°;? 


3.4 ± 1.3 


285.11?^ 


q 1+18.5 

3 - 1 -17.8 


6 < AN < 8 



TABLE I: The best-fit and la confidence level for the velocity dispersion a,, the magnitude (u) and the direction (I, b) of the 
tilted velocity, and the 2a constraints on the number of e-folds of inflation. The "x" means that the value is less than zero. 



Universe" effect, which arises in inflationary models due 
to the perturbations of quantum fields other than the 
inflaton, such as the axion, whose energy density is sub- 
dominant compared to that of the inflaton. There are 
strong constraints on sub-horizon isocurvature perturba- 
tions and thus if they are responsible for the tilt it would 
require them to be much larger on very large scales than 
they are on smaller scales. Although multi-field mod- 
els will not generically produce such a sharp breaking in 
scale invariance or result in isocurvature modes in the 
later Universe, there are double inflation models which 
can achieve this [3] • 

The dipole anisotropy is associated with the isocurva- 
ture fluctuation as follows [3, H| 2 



where 5tp is the fluctuation of the quantum field, and ipo 
is the background field. An inflationary scenario that re- 
sults in isocurvature perturbations in the later Universe, 
such as in Ref. is required. Assume that at the 

onset of inflation, the isocurvature quantum fluctuation 
satisfies 8ip/(po ~ 1 (Constants of 0(1) won't affect the 
results much), thus the Hubble horizon at the onset of 
the inflation is p — 0, I — if,- , and L = e AN H^ 1 , where 
AA~ = N — N m i n is the excess number of inflationary 
c- folds beyond N m i n . The constraints on a "tilted Uni- 
verse" lead to the constraints on the number of e-folds 
of inflation u/c ~ H^ 1 / L ~ e~ AN . Therefore, an ob- 
servable "tilted Universe" velocity requires that inflation 
should last a modest number of e-folds, which should not 
be too long nor too short — if inflation lasts too long, 
such perturbation effects would be washed out, if infla- 
tion is too short, the dipole effect will be too large, mak- 
ing the Universe over-tilted. We show the constraints on 
the number of e-folds in the last column of Table fl] We 
find that the required extra number of e-folds is at least 
6 for all catalogs, and SN, SMAC and COMPOSITE can 
also provide a 2a upper bound on AA~ given their data. 



2 There may be some factors of order unity in this equation for 
different models, but they have negligibly small effect on the 
AiV constraints due to the exponential. The same holds for 
Stp/tpo ~ 1. 



Note that our study provides a similar constraint to the 
Grishchuk-Zel'dovich effect (which can arise from either 
adiabatic or isocurvature perturbations) which requires 
that the extra number of e-folds should be greater than 
about 7 at the 95% confidence interval |22j. Also, corre- 
lations in the CMB multipoles may be used to estimate, 
from the PLANCK data, the tilt with an error bar simi- 
lar to what we obtained here [23j. In the future, there is 
the potential to use the large number of SNe measured 
by the Large Synoptic Survey Telescope (LSST) to con- 
strain u with a standard deviation of about 30 km/s [24j . 
So it may be possible to explore the duration of inflation 
and the pre-inflationary quantum state at quite a precise 
level. The peculiar velocity field, is therefore a powerful 
tool to probe the very early Universe in a manner not 
accessible by CMB observations alone. 

V. CONCLUSION 

In this paper, we developed a model and a statistical 
method to justify whether the apparent bulk flow motion 
of galaxies in our surveys is due to the subtraction of 
the intrinsic dipole on the CMB sky. In the conventional 
bulk flow scenario, the galaxies in the local region (< 150 
/i _1 Mpc) are moving towards direction (I — 287° ± 9°, 
b = 8° ± 6°) with v = 407 ±81 km/s seems in tension 
with the ACDM predictions. However, we point out that 
some fraction of the CMB dipole can be intrinsic due to 
large scale inhomogeneities generated by preinflationary 
isocurvature fluctuations, so that in the CMB rest frame, 
all of the galaxies have streaming velocity towards the 
particular direction, resulting in the tilted Universe. 

We modeled an intrinsic CMB dipole as a tilted ve- 
locity u and developed a statistical tool to constrain its 
magnitude and direction. We found that (1) the mag- 
nitude of the tilted velocity u is around 400 km/s, and 
its direction is close to what was found in previous bulk 
flow studies. For SN, SMAC and COMPOSITE cata- 
logs, u = is excluded at about the 2.5ct level, which 
can explain the apparent flow; (2) there is little corre- 
lation between the tilted velocity u and galaxy's small 
scale velocity and intrinsic dispersion <r* , which confirms 
our assumption. 

Furthermore, assuming that primordial isocurvature 
modes lead to an intrinsic dipole anisotropy, constraints 



6 



on the magnitude of the tilted velocity can result in con- 
straints on the duration of inflation, due to the fact that 
inflation can neither be too long (no dipole effect) nor 
too short (very large dipole effect). Under this assump- 
tion, the constraints on the tilted velocity require that 
inflation lasts at least 6 e-folds longer (at the 95% con- 
fidence interval) than that required to solve the horizon 
problem. 

Finally, we should point out that if there is an intrinsic 
fluctuation, the free electrons in clusters should be able 
to "see" the modulation on the sky, which can be tested 
through the kSZ effect. Therefore, the results from South 
Pole telescope (SPT) [25| and Atacama cosmology tele- 
scope (ACT) [Is H3] may shed some light on the intrinsic 
fluctuations. Such work is in progress. 
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Appendix A: Analytic formulae of correlated 
Window Function } '12 (fc) 

The correlated window function / \ 2 {k) 

/cPk 
— (rrk) (?2-k) x exp (ifck- (n - r 2 )) , 

- 7 . (A1) 

can be calculated analytically, by transforming it into 

harmonic space and using the property of spherical har- 
monics. The final integral should only depend on: (a) the 
angle a between ri and r 2 (therefore ri and r 2 should be 
symmetric) ; (b) r\ and r 2 (amplitude of the vector) ; (c) fc 
(k's amplitude). Therefore, we can specify r\ = (0, 0, 1), 
r 2 = (0, sin a, cos a), where a is the relative angle be- 
tween r\ and r 2 . Therefore, 

A = ri — r 2 = (0, — r 2 sin a, r\ — r 2 cos a). (A2) 

So its direction and amplitude become 

A = — (0, — r 2 sin a, r\ — r 2 cos a), (A3) 



and 



We can set 



A = [r\ + r\ - 2r x r 2 cos a] 5. (A4) 



therefore, 



k • A =— ((— r 2 sin a) sin 9 sin tfi + {r\ — r 2 cos a) cos#) . 

Now we can use spherical harmonic function Yi m (9, (f) to 
decompose the integrand as follows. 

^fi-k^ ^f 2 -k^ = cos6>(sinasin6>sin(/> + cosacos6*) 
IP 



V 15 ^"^i + i 2,-i) 

4 pK I 
+ -J -cos aY 2fi + - cos a, (A7) 



cxp(*fck • (n - r 2 )) = *'(2i + l)ji(kA)Pl$ ■ A), 

(A8) 

in which we can just consider I = 0, 2 two terms in the 
summation. The reason is that the mixing angle be- 
tween k and A just causes the mixing between differ- 
ent m modes in the spherical harmonics in Eq. (IA8|) . so 
the final non I — and I = 2 terms vanish due to the 
orthogonality. Therefore, Eq. (|A8I) becomes 

exp(ifck • (n - r 2 )) = j a (kA) - 5j 2 {kA)P 2 {k • A), 

where ji(kA) is spherical bessel function. 

P 2 (k-A) = i(3(k.A) 2 -I 
3 /2tt 

= -2^2Yyg(r 2 sina) 2 (y 2)2 +y 2 ,_ 2 ) 

+ (( ri -^cosa) 2 - i(r 2 sina) 2 ) r 2 „ 

3 

- -p? (r 2 sin a) (r\ - r 2 cos a) 
A z 

27 



x (y 2 ,i+y 2 ,_i) 



(AI0) 



k =(sin6*cos(/),sin0sin(/), cos 9), (A5) 



Note that there is no zero order term in P 2 (k • A). 
Then we use the orthogonality property of Y\ m 
J (P£Y lm Yfi m , = 5 w 6 mrn , and Y{ m = Y h _ m {-l) m and 
get the final result 

/ia(fc) = ^cosa{j (kA)-2j 2 (kA)) + -^j 2 (kA) ri r 2 sin 2 a. 

(All) 

It is clear that this integration has the three properties we 
listed above and the window function f\ 2 depends only 
on (r±,r2,k,a). This is an independent and simplified but 
equivalent result to Eq.(9.32) in [Hj]. 
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